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Abstract 



This paper studies the performance of antenna array processing in distributed multiple access 
networks without power control. The interference is represented as a Poisson point process. Desired 
^ and interfering signals are subject to both path-loss fading (with an exponent greater than 2) and to 

J> ' independent Rayleigh fading. Using these assumptions, we derive the exact closed form expression for 

(N ; 

the cumulative distribution function of the output signal-to-interference-plus-noise ratio when optimum 
combining is applied. This results in a pertinent measure of the network performance in terms of 
the outage probability, which in turn provides insights into the network capacity gain that could be 

o 

achieved with antenna array processing. We present and discuss examples of applications, as well as 
some numerical results. 
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I. Introduction 

In wireless communications, energy and spectrum resources are scarce and expensive, and 
must consequently be managed efficiently in order to address the growing requirements of 
emerging applications. Intensive work has been dedicated to developing advanced processing 
technologies to improve the spectral efficiency. In decentralized multiple access systems, such 
as ad hoc networks, interference is the major performance inhibitor, which explains why in recent 
years, much research effort has been dedicated to the interference mitigation. Several solutions 
have been proposed involving exploitation of the particular structure of the interference. In fact, 
technologies such as spread spectrum and multiple antennas can be used to design systems with 
a number of degrees of freedom that, when properly exploited, increase system tolerance to 
interference. In the context of spread spectrum, various multiuser receiver schemes have been 
introduced OX The well-known minimum mean square error (MMSE) receiver is the linear 
multiuser receiver that maximizes the signal-to-interference-plus-noise ratio (SINR), and the 
equivalent solution in the context of antenna array processing is known as the optimum combining 
(OC) receiver ||2). Though the study of the reliability of these receivers is an important issue, 
which has been the subject of extensive work S-flU, there is still a lot of work left to do in order 
to understand their limits in various contexts. In fact, the bulk of the results available deal with 
special assumptions concerning interference, such as considering equal power interferers, strong 
interferers or asymptotic cases (infinite number of interferers and antennas). In this paper, we 
address the issue of quantifying the performance gain that could be achieved by employing the 
optimum combining receiver in a decentralized network without using power or access control. 

A. Problem statement and contribution 

As previously mentioned, we are interested in distributed random access networks. A common 
representation of this kind of networks is to consider a random number of interferers that are 
independently and identically distributed over the area. The number of interferers is measured by 
one parameter, namely, their spatial density. Thus, we adopt the so-called stochastic geometric 
model which is widely used in the literature for the access layer as well as for the physical layer 
0-flH. In this model, the link outage is defined as the probability that the expected SINR seen 
by a representative receiver is below a certain threshold. The expectation is taken over the set of 
possible realizations of the network. As defined, the outage probability serves to derive several 
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spatial performance metrics, namely, the mean number of transmissions in outage per unit area, 
the mean throughput per unit area, and the mean distance covered by all transmissions per unit 
area. From these spatial performance metrics, we can study the connectivity, the transport capacity 
and several important issues of the network performance. The tutorial paper IfTOll provides a 
survey of these issues and discusses some applications of this model and this methodology to 
wireless communication. 

To derive the outage probability, we need the interference distribution. In a stochastic geometric 
network, the interference is a spatial shot-noise process with an impulse response having a 
decaying power law form. Even with one antenna and the matched filter receiver, the distribution 
of the interference does not take a closed form (except for some special cases, such as a path 
loss exponent equal to 4) although its characteristic function is well known [fTTTl - [fT4ll . 
Recent works have focused on analyzing the performance of more advanced receivers than the 
conventional one. In [[8]|, lfi~5l and lfi~6l . approximations and bounds on the outage probability 
are given for the successive interference cancellation receiver. In the context of smart antennas, 
ifTTl gives lower and upper bounds on the outage probability for sectorized antenna, maximal 
ratio combining and space time coding techniques. Analyses of the zero-forcing and the partial 
zero-forcing receivers follow in [TT8l . |fi~9ll . The main contribution of the work presented in this 
paper is the exact derivation of the outage probability for the optimum combining receiver with 
an arbitrary number of antennas. 

B. Organization of the paper 

The outline for this paper is as follows. In Section II, we present the system model. The 
derivation of the outage probability, which is the main result obtained in this paper, is presented 
in Section III. Applications and simulations follow in Sections IV and V, and section VI concludes 
the paper. 

II. System Model 

A. Network and channel model 

In a distributed network, a receiver is surrounded by a number of undesired source nodes that 
transmit on the same medium and in the same time slot as its intended transmitter. Usually, 
receiver nodes do not control the number and the positions of these sources of interference; 
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rather, these are determined by the network dynamics and the access strategy in place. In the 
particular case of a random access strategy, a simple way to model the network is to consider 
interfering nodes distributed according to a homogeneous Poisson point process (PPP) Q, [|20ll . 
A PPP relies on a single parameter A, namely, the density of transmitters. In a planar network, it 
is described as follows: In a closed region with area A, the number of transmitters is distributed 
according to the Poisson law with density XA. The positions of these nodes are uniformly 
distributed on the plane. We assume that each node transmits with a single antenna and receives 
with L antennas. The distance between an emitting node and its intended receiver is set to d r . 
Formally, the set of emitting nodes forms a homogeneous PPP II = {Xi 6 ffi 2 , A, i G N}, where 
Xi are the nodes' locations. The channel is modeled by two independent components, the first 
of which represents the path-loss attenuation with an exponent a > 2, and the second, represents 
the channel coefficients, which are independent among antennas and nodes. 

B. Interference expression under antenna array processing 

Since the network is modeled as a homogenous PPP, the interference distribution does not 
depend on the spatial locations of receivers. Thus, in the following we focus, without any loss 
of generality, on a representative receiver placed at the origin. A receiver is subject to a number 
of interfering signals coming from non-desired transmitters. The received signal vector is then: 



where c r and c k , fceN, are the propagation vectors with dimension L that have independently 
and identically distributed, zero-mean, unit variance complex Gaussian entries. The vector n is 
a zero-mean complex Gaussian noise with variance a 2 per dimension. All nodes are supposed 
to use the same transmission power, normalized to unity. 

In statistical antenna array processing, a weight vector is chosen based on the statistics of the 
data received and optimized under a given criterion. Hence, the output of the antenna array 
processor is: 




(1) 



y = w x, 



(2) 
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where w is a complex weight vector. The operator T denotes the transpose conjugate. Several 
criteria could be considered in determining this vector 0, such as a maximization of the desired 
signal-to-noise ratio corresponding to the maximal ratio combining receiver. In systems with 
interference, the aim is to optimize the antenna array output such that the quantity of noise and 
interference is minimized in the resulting signal. In other words, maximizing the SINR: 

/ w T c r c^w \ 
woe = arg max <^ — -— — ■ — > , (3) 

where I L is the L x L identity matrix and Hj is the interference covariance matrix expressed 
as R/ = J2x k en l^fcl ° c fc c fc- The well-known solution corresponds to [0, ETTl : 

w oc = (R/ + ^IlYW (4) 
The resulting SINR is denoted as (3 and expressed as: 

= d~ a c^R _1 c n (5) 
where R is the interference plus noise covariance matrix expressed as R = R/ + <t 2 Il. 

III. Outage probability derivation 

The configuration we consider has two random parameters, namely, the locations of interferers 
and the channel coefficients. Thus, we have a complex probabilistic system, which we propose to 
break down into two levels. First, we generate the SINR expression conditioned on a realization 
of the network, after which we average over this random quantity. Consider an arbitrary planar 
region D with finite radius d and arbitrary nodes positions, and use iV to denote the number of 
nodes in place. The problem now is to find the distribution (or the eigenvalue distribution) of 
random quadratic matrices having the form (CPC T + <j 2 \l). Each column of C represents a 
channel vector of an interferes Thus, the random matrix C has column vectors independently 
and identically distributed as multivariate normals with zero-mean and unit- variance vectors. The 
matrix P = diag\\Xi\~ a , |X 2 | _Q , • • •] is a diagonal matrix with real elements corresponding to 
the set of received powers. 
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A. Related results 

The problem stated in the previous paragraph is a classic problem in the probability and 
wireless communication literature. We can classify the most pertinent results to the context of 
this work into two categories. The first category concerns the work of Silverstein and Bai ll22ll . 
who consider an asymptotic regime (the number of nodes and the number of antennas tend to 
infinity at a fixed rate). They established that the eigenvalues of the considered class of random 
matrices converge to a deterministic limiting distribution. This result is first used in ll23l and 
ll24l to derive a closed form expression of the asymptotic mean and variance of the SINR at the 
output of the MMSE receiver, in the case of deterministic received powers. In [25J, the authors 
apply these asymptotic techniques to the particular case of emitting nodes uniformly distributed 
on an infinite plane, and provide the SINR mean and variance. Following this methodology, ll26l 
gives an approximation of the SINR distribution based on the assumption that the latter is a 
Gamma distribution. In our case, despite the large number of users (the network area is wide), 
these results cannot be applied. In fact, the spatial separations of users ensures that the global 
interference is only influenced by a small number of them [H, D3L lfl6ll . In addition, the goal 
of our work is to show how using a small number of receive antennas can improve the SINR. 
The second category was initiated by Khatri ll27ll . who derived the distribution of the matrix 
CPC T . From this, ll28l provides the SIR distribution when the noise is ignored: 



where Hq is a hypergeometric function of matrix arguments, Z = diag[(l + qfi)~ x , It,-i], q is 
a particular constant and T(-) is the gamma function. 

Even if we consider that our system is interference-limited, expression © is not simple as it 
contains the hypergeometric function with matrix arguments. It therefore appears to be difficult 
to derive the SINR distribution when P is a random matrix. In [|29ll and [30], the authors 
extended the previous result by including the noise term and derived a simpler expression. 
The derivation consists in expressing the hypergeometric function as a polynomial ratio. The 
cumulative distribution function (CDF) of the SINR is given as: 



r(iV + l) (3 L - 1 q N+l 



P]- 1 ^ r } (N+V,I N -qP-\Z 



(6) 



T(L)T(N+1-L) (l + q/3) N + l 



F f) (P\N,X 1 ..-X N ) = l 



(7) 



exp (a 2 /3d?) nf=i(l + \Xj\- a Pd?) ' 
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B. Outage probability expression 

First, let us consider the case of a single receive antenna that corresponds to a Rayleigh 
fading channel. In this case, as established by 0, the complementary CDF of the SINR is 
the characteristic function of the sum of interference and noise, which is a perfect match for 
expression ©. In fact, rijli(l + \Xj\~ a l)~ l ls me characteristic function expression of a 
weighted sum (with weights equal to |X, |~ Q ) of N independently and exponentially distributed 
random variables. Since the characteristic function of a Poisson class of interferers is well known 
ifTT fl -lfH l - [|20l . OTTl . |[32ll . the derivation for L = 1 is trivial. For an arbitrary number of antennas, 
taking the expectation over the possible realizations of the network and denoting (3d" by 7, we 
establish the following theorem. 

Theorem 1: Using the MMSE receiver, the outage probability in a Poisson field of interferers 
and Rayleigh fading channel is: 

F(7 , A) = ! _ (A^1±£V exp( _ AV /« _ ^ 

(8) 

where A = 7r2/«r(2/«)r(l - 2/a). 

Proof: The coefficients cij in © can be derived simply. In fact, extending the exponential 
in the denominator of © and putting Pj = \Xj\~ a we get: 

ex P ( ( x 2 7 ) + ^7) = £ E ^( p i> ■ ' ' • (9) 

j=l fc=0 ' i=0 

where &i(Pi, ■ ■ • , Pjv) = Si<i 1 <-<i i <jv ^h^h ' ' ' ^k- Then, ai are expressed as follows: 

i / 2\i— k 

ffli = E7Tw 6t ( Pl ''"' Pw ) z = l---L-l, (10) 

fc=0 ^ 



From © and CLS): 

~t 1 «1 r 



F 7 ( 7 ,A) = 1 - exp(-a i 7 )E 



N 



E E £V^- -» 

i=0 fc=0 v y 



nf=i(i+^-7) 



(id 
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where E x denotes the expectation with respect to the random variable x. The coefficients 



h(Pi 



P N ). k 



1---L — 1, are composed of a sum of products of combinations of k 



elements from the set {Pj,j = 1 ■ • ■ iV}. Given that the locations of the nodes are independently 
and identically distributed, the last expression simplifies to yield: 



E 



Pi,-,Pn 



b k (P u --- ,P N )^ 

nL (i + iVr) 



Pl---Pkl k 



1 



C N E Pl 



Pil 



E 



Pi 



[1 + Pij) 



N-k 



(12) 



where is the number of combinations of size k from a set of N elements. To compute 



E 



Pi 



Pi 



recall that conditioned on N, in a finite region D 6 M. 2 , the locations of nodes 



a+pu) 

are uniformly distributed. Consequently 



E 



Pi 



\ P1 




L(i+p l7 ). 





Similarly, we have: 



E 



Pi 



1 


= 1 f 


L(l + Pi7). 


■nd 2 J D 



;i + \x\- a j)- l dx. 



(13) 



(14) 



Since the number of nodes is distributed according to the Poisson law with mean Xnd 2 , and 
considering ([13]) and ([I4i the CDF of the SINR yields: 



F 7 ( 7 ,A) 



oo L— 1 min(i,N) 

l-exp(-a 2 7 )EE £ 

N=0 i=0 k=0 
N-k 



nd 2 



dX 



D 



1 + |x|- a 7/ 

L— 1 i oo 



k\(N-k)\(i - fc)! 
(\nd 2 ^ N 



nd 2 



D 





X 


\~ a l 


1 + 1 


X 


-dry 



N\ 



exp(—\ird 



i-«p(-^7)EEE 



(^ 2 7) 



1 



(N — k)\ 



A 



i=0 fe=0 iV=fc 

1 + |X|- Q 7 



fc!(i - A;)! 



A 



D 







~ Q 7 


1 + 1 







dx 



D 



N-k 



exp(—Xnd 2 



Using the fact that nd 2 = J D ldX, we get: 



L-l 



^ 7 ( 7! A) = l-exp(-a 2 7 )^^ 



2_\i-k 



i=0 k=0 



kl(i - k)\ 



(^7) 



!X| "5-<1X) ,x P (A / ( 



1 + |X|- Q 7 



-l)dX). 
(16) 
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Finally, letting the area of D go to infinity and evaluating the integrals in (fT6l) (the complete 
evaluation of these integrals is provided in the Appendix) we get: 

L-1 i 1 

F 7 ( 7 ,A) = 1 - ex P Hx 2 7 ) kUl _ k )M 2 ^ l ~ k (7 2/a AA) fc exp(- 7 2 /°AA) 

i=0 fc=0 

^ (AA 7 2 /" + a 2 iY , 2/a 2 . 
= 1~2^ exp(-AA 7 ' -(T7)- 

8=0 

(17) 
■ 

IV. Discussion and application 

A. Noise-limited regime and interference-limited regime 

The expression ([8]) derived for the outage probability has a simple form. From it, we can easily 
see the effects of various parameters on the network performance. Moreover, it clearly presents 
the trade-off between noise and interference cancellation, and thus provides some insight into 
the intuition behind the result. 

1 ) Noise-limited regime: In systems where the density of users is negligible, the SNR cumu- 
lative distribution function becomes: 

F 7 ( 7 ) = 1 - ^r^xp(-a 2 7 ). (18) 

i=0 

As expected, (fl"8l is the classic expression of a Chi-square cumulative distribution function. 
The latter provides the SNR distribution when maximal ratio combining is performed, that is 
equivalent to optimum combining in a complex Gaussian noise environment. 

2) Interference-limited regime: Where the noise is negligible, the CDF is: 

F 7 ( 7 , A) = 1 - I ex P (-AA 7 2/a ). (19) 

i=0 

Equation (fT9l) corresponds, up to the factor A, to the probability that the L th largest received 
power is below the threshold 7 [33 J. This is not surprising given some properties of the interfer- 
ence. Since the path loss exponent is greater than 2, the interference distribution is heavy tailed 
ll34l . Therefore, there is a large dispersion between the received powers. Let us order the received 
powers and denote them as |X(x)|~ a 3> |X( 2 )|~ Q 3> |X( 3 )| -a 3> • • • . With probability one, the 
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random matrix CPC T is of order L, and its eigenvalues could be approximated as |X(j)|~ a , 
i — 1 • • • L ll35ll . Thus the maximal eigenvalue of the inverse of CPC T is approximately \Xr L \ \ a . 
The performance is then primarily affected by the L th strongest received power, and so, the outage 
is almost due to the L th strongest interferer, and corresponds to the event \X(i,\\ < 7 1//a . In other 
words, to achieve successful reception, the L th strongest interferer must be outside the region 
of radius 7 1 /". This fact matches expression (TT9l) up to the scalar factor A. 

3) SIR mean and variance: From (fl9l ), we can establish in a straightforward manner that in 
an interference-limited regime, the mean and the variance of the SIR are, respectively: 

E[SIR] ~ (L-l)! (20) 



Var[SIR] 



T(L + a) 



r(L + a/2) 
(L-l)! 



^ d~ 2a 



r 



\ a A c 



(21) 



(L-l)! 

Thus, the MMSE receiver provides an antenna array gain equal to T(L + a/2)/(L — 1)!. It 
should be noted that expression (fl9l) is consistent with the result previously established in [|25l . 
In fact, when the number of antennas is sufficiently high, the mean SIR is approximately equal 
to L a l 2 xa f.2^ a/2 ■ This latter relation is also derived for the asymptotic regime, i.e., very large L, 

m ma. 



B. Application: single-hop throughput capacity 

The interference model considered serves to capture several network classes under some 
additional assumptions. In the following, we focus on a single-hop ad hoc network with an 
ALOHA access protocol. Thus, the density of interferers represents the spatial rate at which 
transmissions occur, i.e., the contention density. In this context, the outage probability represents 
the spatial average of the density of communications that fail to be established at a given range 
d r in almost every given realization of the network Q. Equally, the mean number of successful 
transmissions or the throughput per unit area is T = A(l — F(j,X)). The simplicity of our 
analytical result allows the direct optimization of the contention density as: 

A m ax = arg max A(l - F(j, A)). (22) 

0<A<oo 

Corollary I: In an interference-limited regime, the optimum contention density is: 

A - 9 ^ (2^ 

AmaX ~ /\ry2/a ■ 
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The parameter g(L) is only a function of the number of receive antennas and corresponds to 
the positive root of the following polynomial: 

L ~ l 4-i 4-L 

Moreover, this parameter satisfies: 

\ < 9(L) < L (25) 

The equivalent spatial throughput per unit area is: 

g L+1 (L) 

Tmax = (L-i)!A 7 2 /° ex P("^( L ))- W 
Corollary \T\ indicates that the MMSE receiver provides a linear scaling of the optimum 
contention density with the number of receive antennas unlike to the MRC and the zero-forcing 
receivers, whose scaling laws are L 2 / a and L 1_2 / Q , respectively ifTTl . [fT8l . The linear scaling 
law achieved by the MMSE receiver is predicted in |[T9l . In fact, the authors show that linear 
scaling is possible with a partial zero-forcing receiver, which is suboptimal as compared to the 
MMSE receiver. Through corollary [T] we confirm this prediction, and we show that the scaling 
law is exactly linear. 

V. Simulation Results 

All simulations are carried out with the following parameters: path-loss coefficient a = 3.5, 
SINR threshold (3 = 3dB and distance between transmitter and receiver d r = 10m. Figure \T\ 
shows the simulation and the analytical results of the outage probability as a function of the 
density A where 1, 2, 3 and 4 antennas are used. The Monte Carlo simulation and the analytical 
curves are very close, with the gap between them arising from the fact that the analysis performed 
on the previous section concerns infinite networks. In our simulations, the area is finite, and 
depends on the density A. It is chosen such that we have 100 emitters in the area, on average. 
In 11361 . the author provides a complete analysis of the error on the interference estimation when 
a finite network is considered rather than an infinite one, and establishes that this error depends 
on the exponent a. 

Figure [2] provides a performance comparison between the optimal combining receiver and three 
antenna array processing techniques, namely, maximal ratio combining, zero-forcing and partial 
zero-forcing receivers. It is clear that the optimum combining receiver outperforms the other 
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techniques. This is due to the fact that maximal ratio combining deals only with the fading effects 
of the desired signal. Zero-forcing receiver uses all additional degrees of freedom provided by 
the antennas to cancel strong interferers, while the partial zero-forcing receiver uses some of 
the antennas for interference cancellation and provides diversity with the remaining ones; they 
are nevertheless still suboptimal as compared to the MMSE receiver, which provides the best 
trade-off between interference cancellation and spatial diversity. 

Figure [3] presents the throughput density improvement with the number of antennas and figure 
H] shows the linear scaling of the optimum contention density. 



This paper derived the exact outage probability of the optimum combining receiver in the 
presence of noise and Poisson field of interferers. The framework developed is nonspecific and 
can be generalized in a straightforward manner to any fading having a power distribution of the 
form J2k xke ~ kx \ such as the Nakagami fading. The result obtained has a simple closed form, and 
provides an understanding of the performance of the MMSE receiver, in addition to allowing a 
comparison with other antenna array processing methods, such as maximal ratio combining and 
zero-forcing receivers. The assumption made on the interference is pertinent to many network 
classes. In addition, the result can be used to study the performance of the optimum combining 
receiver at the access and at the network layers. Finally, the simulations provided in the paper 
demonstrate that experimental and theoretical results match perfectly. 

Appendix A 
Evaluation of the integrals in (TT61) 

The integrals in (fT6l) are evaluated as follows: 

/ ((1 + (AT^)- 1 - l)dX = / / ((l+r- a 7) -1 -l)rd0dr 



VI. Conclusion 





(27) 



Putting y = r a 7, the last expression becomes: 
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((1 + ixi"^)- 1 -l)dX 



7 



« J o 

2/a 



1 + y 



-2tt- r(2/a)r(l - 2/or) 

a 

- 7 2 / Q A. 



(28) 



On the other hand, we have 





X 


\~ a l 


1 + 1 


X 


-dry 



dx 



((1 + IXI"^)- 1 - l)dar 



(29) 



Appendix B 
Proof of corollary 1 

By Descartes' rule of signs, the polynomial in expression (l24l) has at most one real positive 
root. The value of this polynomial at L is negatively signed: 

L L 



Q(L) 



i=0 



(L-l)! 



i=0 
L-l 

= E 

i=0 
< 



Z7 



(30) 



On the other hand, the value of the polynomial at L/2 is lower bounded as follows: 

L L 



Q(L/2) 



L-l 

Y — 



i=0 
L-l 



> 



> 



2H\ 2 L (L-l)\ 
L L 



i=0 

3 
2 



2* 2 L (L-1)! 



1 



2 L 2 L (L-1)! (31) 
The right hand side of inequality (|3TT) is positive for all values of L. Thus, the considered 
polynomial has one positive root that lies in the interval [L/2, L\. 
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Fig. 1: Outage probability: simulation (dashed lines) and theoretical (solid lines) results with 
a 2 = —50dB / antenna (all nodes use unit transmission power.) 
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Fig. 2: Outage probability: Maximal ratio combining, zero forcing, partial zero forcing (with the 
optimum number of cancelled interferers established in |fl9ll ) and optimum combining with L=3 
and a 2 = 0. 
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Fig. 3: Optimum combining: Throughput spatial density with L = 

—57dB / antenna. 
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Fig. 4: Optimum combining: Optimum contention density as a function of the antenna array 
dimension, with a 2 = and A7 2 / a normalized to 1. 



